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ABSTRACT 

We study future constraints on dark energy parameters determined from several com- 
binations of cosmic microwave background experiments, supernova data, and cosmic 
shear surveys with and without tomography. In this analysis, we look in particular for 
combinations of experiments that will bring the uncertainties to a level of precision 
tight enough (a few percent) to answer decisively some of the dark energy questions. 
In view of the parameterization dependence problems, we probe the dark energy using 
two variants of its equation of state u'(z), and its energy density pde{z). For the latter, 
we model pde{z) as a continuous function interpolated using dimensionless parameters 
£i{zi) = pde{zi) / pde{0) . We consider a large set of 13 cosmological and systematic 
parameters, and assume reasonable priors on the lensing and supernova systematics. 
For CMB, we consider future constraints from 8 years of data from WMAP, one year 
of data from Planck, and one year of data from the Atacama Cosmology Telescope 
(ACT). We use two sets of 2000 supernovae with Zmax = 0.8 and 1.5 respectively, 
and consider various cosmic shear reference surveys: a wide ground-based like survey, 
covering 70% of the sky, and with successively 2 and 5 tomographic bins; a deep space- 
based like survey with 10 tomographic bins and various sky coverages. The one sigma 
constraints found are {a{wo) = 0.086, (j{wi) — 0.069}, {cr(u'o) — 0.088, a{wa) ~ 0.11}, 
and {ct{£i) = 0.029, CT(f2) = 0.065} from Planck, supernovae and the ground-based 
like lensing survey with 2 bins. When 5 bins are used within the same combination the 
constraints reduce to {(t(wo) — 0.04, a{wi) — 0.034}, {ct(wo) — 0.041, (T(u'a) — 0.056}, 
and {(j{£i) = 0.012, ct(£^2) — 0.049}. Finally, when the deep lensing survey with 
10% coverage of the sky and 10 tomographic bins is used along with Planck and the 
deep supernovae survey, the constraints reduce to {a{wo) = 0.032, (t(wi) = 0.027}, 
{a{wo) = 0.033, cr(wa) = 0.04}, and {^(fi) = 0.01, ^(fz) = 0.04}. Other coverages of 
the sky and other combinations of experiments are explored as well. Although some 
worries remain about other systematics, our study shows that after the combination 
of the three probes, lensing tomography with many redshift bins and large coverages 
of the sky has the potential to add key improvements to the dark energy parameter 
constraints. However, the requirement for very ambitious and sophisticated surveys in 
order to achieve some of these constraints or to improve them suggests the need for 
new tests to probe the nature of dark energy in addition to constraining its equation 
of state. 

Key words: cosmology: theory - dark energy - gravitational lensing - large-scale 
structure of universe 



1 INTRODUCTION 

One of the most important and challenging questions in cosmology and particle physics is to understand the nature of the 
dark energy that is driving the observed cosmic acceleration, see e.g. ( [Weinberg (1989)1 [CarroU et al. (1992)] [Turner (2000)| 
|Sahni fc Starobinsky (2000) [ [Padmanabhan (2003) [ [Ishak (2005)[|. An important approach to this problem is to constrain 
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the properties of dark energy using cosmological probes. This would provide measurements that would allow one to test 
various competing models of dark energy. However, due to a high degeneracy within a narrow range of the parameter space, 
constraining conclusively dynamical dark energy models is going to be a difficult goal to achieve, and much effort and strategy 
will be needed. Ultimately, a combination of powerful cosmological probes and tests will be necessary. 

In this paper, we study how dark energy parameters are constrained from different combinations of cosmic microwave 
background (CMB) experiments, supernovae of type la (SNe la) data, and weak lensing surveys (WL) with and without 
tomography. In particular, we look for combinations of experiments that will be able to constrain these parameters well 
enough to settle decisively some of the dark energy questions, say to a few percent. When CMB measurements constrained 
the total energy density to D,t = 1.02 ± 0.02 to a one sigma level jSpergel et al. (2003) | [Rennet et al. (2003)) , it became 
generally more accepted that spatial curvature is negligible. Thus, an uncertainty of a few percent on dark energy parameters 
could be set as a reasonable goal. Of course, one should bear in mind that it will remain always possible to construct dynamical 
dark energy models that are indistinguishable from a cosmological constant within these limits, and therefore, one needs to 
resort to cosmological tests beyond the equation of state measurements. A better scenario providing a decisive answer would 
be one in which one could show that dark energy is clearly not a cosmological constant. 

It is certainly wise to probe the nature of dark energy using gradual steps. However, in both the scenarios mentioned 
above, one should keep in mind that the results and conclusions obtained from an analysis where the equation of state is 
assumed constant are subject to changes if the equation of state is allowed to vary with the redshift. In this paper, we consider 
dark energy with a varying equation of state. 

We chose the combination CMB+SN la+WL as various studies have already shown that supernovae type la constitute a 



powerful probe of dark energy via distance-redshift measurements, see for example ( 


Riess et al (1998)| 


Garnavich et al. (1998)| 


IFilippenko & Riess (1998)||Perlmutter, et al. (1999)| 


IPerlmutter S., et al. (1997)| 


Riess et al (2000)| 


Riess et al (2001)|l 


JTonry et al. (2003) | |Knop et al. (2003) | 


IBarris et al„ (2004)| |Riess et al. (2004) 


. Also, several parameter forecast studies 



have shown that combining constraints from weak gravitational lensing with constraints from the CMB is a powerful combina- 



tion to constrain dark energy; see, e.g. l|Hu (2001)||Huterer (2002)||Huterer & Turner (2001)| 


Benabed & Van Waerbeke (2003) } 


Abazajian & Dodelson (2003)||Refregier et al. (2003)| [Heavens (2003) | 


[Simon et al (2003) 


|Jain & Taylor (2003) |l 


Bernstein & Jain (2004)||Song & Knox (2004)|l. Importantly, weak lensing measurements are sensitive to both the effect of 



dark energy on the expansion history and its effect on the growth factor of large-scale structure. Furthermore, in addition 
to tightening the constraints, using independent probes will allow one to test the systematic errors of each probe, which are 
serious limiting factors in these studies. For each of these probes, much data will be available in the near and far future. 
For WL, there are many ongoing, planned and proposed surveys, such as the Deep Lens Survey (http://dls.bell-labs.com/) 
jWittman et al.(2002)t ; the NOAO Deep Survey {http:// www.noao.edu/noao/noaodeep/ ); the Canada-France-Hawaii Tele- 
scope (CFHT) Legacy Survey (http://www.cfht.hawaii.edu/Science/CFHLS/) jMellier et al.(2001)^ ; the Panoramic Sur- 
vey Telescope and Rapid Response System { http://pan-starrs.ifa.hawaii.edu/ \; the Supernova Acceleration Probe {SNAP; 
^tp://snap.lbl.gov/\ i |Rhodes et al.(2003)||Massey et allW OS) Rcfr egier et al.(2003)| l; and the Large Synoptic Survey Tele- 
scope (LSST; http://www.lsst.org/lsst_home.html\ | ^ysoi^(2002)(). Similarly, there are many ongoing, planned and proposed 
SNe la surveys, such as the Supernova Legacy Survey jPain et al. (2002) | [Prichet (2004)11 (SNLS); The Nearby Super- 



nova Factory (SNfactory) JWood-Vasey et al. (2004)] |; th e ESSENCE project jSmith et al. (2002) | [Garnavich et al (2002)| 
[Kirshner et al. (2003)t; Sloan Digital Sky Survey (SDSS) JMadgwick et al. (2003)t; Th e Carnegie Supernova Project 
jFYeedmeonst al. (20 04)1 1; and the Dark Energy Camera Project~ ]DeCMn^r£^2004)[ |. We should mention that there are 
other noteworthy cosmological tools for probing dark energy that we did not consider in this study, notably clusters of 
galaxies (see for example ( Mohr (2004)[ [Wang et al.(2003)1 l and references therein), Lyman-alpha forests (see for example 



jMandelbaum et al (2003)"] Se^a^^^^^20M^| l), and baryonic oscillations (see for example ( [Eisenstein (2003)|[Seo fc Eisenstein (2003)[ 
Blake fc Glazebrook (2003)[[Linder (2003)| l). 



For the CMB, we consider future constraints from 8 years of data from the Wilkinson Microwave Anisotropy Probe 
(WMAP-8) | [Bennet et al. (2003)[|Spergel et al. (2003)^ , 1 year of data from the Planck satellite (PLANCKl), and 1 year of 
data from the Atacama Cosmology Telescope (ACT), see e.g. ( [Kosowsky (2003)1 l. We use two sets of 2000 supernovae with 
Zmax=0.8 and Zmax=1.5 respectively, and consider two types of cosmic shear surveys: a ground-based like survey covering 
70% of the sky with source galaxy redshift distribution having a median redshift Zmed ~ 1, and a space-based like deep survey 
covering successively 1%, 10% and 70% of the sky with Zmed ~ 1.5. 

We take into account in our analysis systematic limits for the supernovae by adding a systematic uncertainty floor 
in quadrature following ( [Kim et al. (2003)] l. We also include for weak lensing the effect of the redshift bias and the shear 

calibration bias by adding and marginalizing over the corresponding parameters as in jlshak et al. (2004) 

The constraints on the dark energy equation of state are parameterization dependent; see, e.g. jWang fc Tegmark (2004)] 
[Upa£lhye_et_aL_(^^ Also, there is a smearing effect due to double integration involved when using the equation of state 



jMaor et al. (2002)[ [Maor et al. (200T) 



In order to partly avoid this, one could probe directly the variations in the dark 



energy density using the data; see, e.g. i Wang fc Mukherjee (2004)] [Wang fc Freese (2004)[ |. However, it has been argued that 



the equation of state is closer to the physics, as it also contains information on the pressure, and, trying to probe the equation 
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of state from the density leads to instability and bias jLinder (2004)| l. Therefore, we choose in this analysis to use both and 
parameterize the dark energy using its density pde{z) as well as two different parameterizations of its equation of state w{z). 



2 COSMOLOGICAL MODEL AND DARK ENERGY PARAMETERIZATION 
2.1 Model 

A total set of 13 parameters is considered as follows. For constraints from WL, we use: Q,rnh^ , the physical matter density; 
Qa,wo and wi (or Wa), respectively the fraction of the critical density in a dark energy component and its equation of state 
parameters (see Sec. 12. 21 (alternatively, we use the dark energy density parameters £i = Pde{zi) / Pde{0) with i = 1,2 (see section 
12.21 '): we use ns(fco = 0.05/i/Mpc) and a^, the spectral index and running of the primordial scalar power spectrum at kg; ag", 
the amplitude of linear fluctuations. In order to parameterize some systematics, we include as a parameter Zp, the characteristic 
redshift of source galaxies (see Eq. I13II , as well as and ("r , the calibration parameters as defined in Ref . I Ishak et al. (2004)| 



which determine the absolute calibration error on the lensing power spectrum, ( [Hirata fc Seljak (2003) l), and the relative 



calibration between tomography bins (see Sec. 13. When we combine this with the CMB, we include Q,i,h^, the physical 
baryon density; r, the optical depth to reionization; T/S the scalar-tensor fluctuation ratio. We assume a spatially flat Universe 
with flm + SIa — 1. This fixes flm and Ho as functions of the basic parameters. We do not include massive neutrinos, or 
primordial isocurvature perturbations. For the supernova analysis, we use Qa, wq, ifi, (or £i, £2) and treat the magnitude 
parameter as a nuisance parameter. We use the fiducial model (e.g. Ref. ( [Spergel et al. (2003)} and add wo and wx): 
Q.bh^ = 0.0224, n„i/i2 = 0.135, = 0.73, wq = -1.0, = 0.0, = 0.93, as = -0.031, as = 0.84, r = 0.17, T/S = 0.2, 
Zp = 0.76, 1.12, Cs = 0.0, and Cr = 0.0. 

2.2 Dark energy parameterization 

As mentioned earlier, constraining the dark energy using its equation of state is known to be parameterization dependent, e.g. 
Wang & Tegmark (2004) JlJpaxlhye^^al. (2004)) , and also suffers from a smearing due to the double integration involved 



Maor et al. (2002) lMaore^d^2o'oi)| |. Alternatively, one can probe directly the variations of the dark energy density as a 



function of redshift pde[z). On the other hand, it has been argued that the equation of state contains information on both the 
density and pressure of the dark energy and using the density to probe the equation of state may lead to instability and bias 
jLinder (2004)| |. We chose to study the constraints on dark energy using both approaches. 

2.2.1 The equation of state 

There are several parameterizations of the dark energy equation of state that have been used to study currently available data 
or to do parameter constraint projections. Discussions of the advantages and drawbacks of some of these parameterizations 
can be found in jWang &: Tegmark (2004")] [Upadhye et al. (2004)| . We used here the following two parameterizations, which 
have no divergence at very large redshift: 

a) (100, wi) 

Here w\ represents the redshift derivative of w^z) in the recent past as follows (see, e.g. ( [Upadhye et al. (2004)| |) 

I If + wiz if 2 < 1 , , 

vj(z) ~ { ^ (1) 

^ ^ 1 Wo + wi if 2 > 1. 

The evolution of dark energy density with redshift is given by pde{z) — pde{0)£{z) where 

. / (i + 2)3(i+™o~™i) e^^i- if2<l, 

I (1 + 2)^(l + "'0+™l) e3™l(l-21n2) if^>]^^ (.^^ 

b) {W0,Wa) 

Here the equation of state is parameterized as jChevalier et al. (2001)| [Linder (2003)| l 

w{a) = wo + Wa ^ = wo + Wa{l — «) (3) 

r ~r z 

where a is the scale factor. The dark energy density evolves with £{a) now given by 
£{a) = a-3(i+»o+»o)g-3»„(i-a)^ 
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Table 1. CMB experiment specifications for Planck and ACT. The parameters used for WMAP 8 years are based on the projection of 
the one year operation and are described in section 5. 
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Figure 1. Convergence auto power spectra 
bottom to the top, the curves correspond to 
errors averaged over bands in I. 



ior the 10 tomography bins. All the parameters are fixed at their fiducial values. From the 
the redshift intervals [0.0, 0.3], [2.7, 3.0]. For each curve, we display the sample variance 



2.2.2 The density parameters: Ex = -77^, £-'2 = ^f^W 

Following jWang fc Mukherjee (2004)||Wang fc Freese (2004)^ , we parameterize £-{z) = ^^^l^j as a continuous function inter- 
polated between today and its amplitude parameters £1 and £2 corresponding respectively to z=0.5 and 1.0, and remaining 
constant at higher redshifts. We use a polynomial interpolation as in | |Wang fc Mukherjee (2004)||Wang fc Freese (2004)| l so 

z 

£(z) = 1 + (4£:i -£2- 3) + '1{£2 - 2£x + 1)^—. (5) 

Zmax Zj^fj^j. 

where the parameters £1, and £2 will be constrained from the data. As suggested in ( |Wang fc Mukherjee (2004)| | and 
jWang fc Freese (2004)] |, we could use more density parameters than two as much more data will be available in the future, 
but we chose to use only two parameters in order to keep the number of parameters equal to the equation of state case and 
be able to make a fair comparison of the results. Departures of the density parameters from unity will be an indication of a 
redshift evolution of the dark energy density and will rule out a cosmological constant. 



3 PROBING DARK ENERGY WITH COSMIC SHEAR 



Weak lensing is a very promising tool for an era of precision cosmology. Already, several studies used currently available cos- 
mic shear data to constrain various cosmological parameters ( |Contaldi et al.(2003)||Van Waerbeke et al.(2002) et al.. (2002) | 
[Wang et al.(2003)]|Jarvi s et al.(2003)||Massey et al.(2004) |. Using statisti cal inference theory, many other studies showed the 
promise of this probelHu fc Tegmark (1999)||Hu (2001)| Huterer (2002)| [Abazajian fc Dodelson (2003)|l 



Benabed & Van Wacrbcko (2003) Takada fc Jain (2004) [Takada fc White (2003) | [Heavens (2003) ||Jain fc Taylor (2003)|l 



Bernstein fc Jain (2004) | [Ishak et al. (2004) | |Sinior^^l(20'03)| |. In particular, weak lensing was shown to constrain signif- 



icantly the dark energy parameters. The advantage of weak lensing is that it is sensitive to the effect of dark energy on 
the expansion history and its effect on the growth factor of large-scale structure. Another advantage of weak lensing is that 
it allows one to construct new tests or techniques to probe cosmology. These include redshift bin tomography | |Hu (1999)| 
|Hu (2002)| |, cross-correlation cosmography jBernstein fc Jain (2004) J , and the use of higher order statistics such as the bis- 
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WMAP8+SN 

WMAPS+SN+WL 
WMAP8+WL 
SN+WL 




Figure 2. The la confidence two-dimensional regions (Ax^ = 1) for the (wq, w\) parameters. The plots contrast the constraints obtained 
from different combinations of 8 years of data from WMAP, 2000 SNe la with Zmax = 0.8, and a ground-based like WL reference survey 
with fsky = 0.7, n = 30gal/arcmin^, (7f„t)^''^ ~ 0.4, and a median redshift z^ed ~ 1- 



pectrum; see, e.g. l |Takada fc Jain (2004)| l. We explore in this analysis the constraints obtained from the convergence power 
spectrum and multiple bin tomography. 



3.1 Convergence power spectrum 

Light rays traveling to us from background galaxies get deflected by mass fluctuations in large scale structures. This results 
in distortions of the sizes and shapes of these galaxies that can be described by the transformation matrix 



del 



K — 7l 

72 



72 

K + 7i 



(6) 



where 6s is the angular position in the source plane; is the angular position in the image plane; k is the convergence and 
describes the magnification of the size of the image; 71 and 72 are the components of the complex shear and describe the 
distortion of the shape of the image. In the weak gravitational lensing limit, |fi;|, I7I <^ 1. 

The convergence is given by a weighted projection of the matter energy density fluctuations 5 = Sp/p along the line of 
sight between the source and the observer. 



<0) 



W{x)S{x,xO)dx 



where x is the radial comoving coordinate and Xh the comoving coordinate at the horizon. 

The convergence scalar field can be decomposed into multipole moments of the spherical harmonics as 

Lm 

where 

1^1^ = J deK{e,x)Yr*ie). 

The convergence power spectrum PI' is then defined by 

{K.lmK,l'm') ~ &U' 5mm' Pi 



(7) 



(8) 



(9) 



(10) 



and we will use it as our weak lensing statistic. In the Limber approximation, it is given by | |Kaiser (1992)|pain fc Seljak (1997)| 
[Kaiser (1998)1 1: 



XH „2 



P-iD 



I 



X dx- 



(11) 



a^ix) Vsinif(x)' 

where P^d is the 3D nonlinear power spectrum of the matter density fluctuation, 5; a{x) is the scale factor; and sin^ x = 
A'-i/2sin(A'i/2^) is the comoving angular diameter distance to x (for the spatially flat universe used in this analysis, this 
reduces to x)- The weighting function g{x) is the source-averaged distance ratio given by 



(12) 
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Figure 3. The la confidence two-dimensional regions (Ax^ = 1) for the {wo,wi) parameters. The plots contrast the constraints 
obtained from different combinations of 1 year of data from Planck, 2000 SNe la with Zmax = 0.8, and a ground-based like WL reference 
survey with fsky = 0.7, n = 30gal/arcmin^, (7f„t)^''^ ~ 0.4, and a median redshift z^^^ 1. 



where n{xiz)) is the source redshift distribution normahzed by f dz n(z) = 1. For n{z), we use the distribution | |Wittman et al.(2000)^ : 

"W^T^e-^/^", (13) 



which peaks at Zp — 2zo. For cosmic shear calculations, we integrate numerically the growth factor using jLinder fc Jenkins (2003)} 



G' 



7 3 w{a) 
2 ^ 2l+X(a) 



G' 3 1 - w(a) G 







(14) 



a 2 l+X{a) a2 
where G — D/a is the normalized growth factor with D — 5(a) / 5{ai) , 

^ (l-an)a,«£:(a)' 

and £(a) is as given in Eq. |1] We use the mapping procedure halofit jSmith et a. (2003) | l to calculate the non-linear power 
spectrum. We show in Fig.Qconvergence power spectra for the 10 tomographic bins. We also show the sample variance errors 
averaged over bands in I. 



(15) 



3.2 Weak lensing tomography 

The separation of source galaxies into tomographic bins improves significantly the constraints on cosmological parameters, 
and particularly those of dark energy because tomography probes the growth of structure. The constraints obtained within 
different redshift bins are complementary and add up to reduce the final uncertainties. We explore here three tomography 
studies using two different types of cosmic shear surveys. For the first survey (ground-based like), the source galaxy redshift 
distribution has a median redshift Zmedian ~ 1.0, and we assume that the photometric redshift knowledge will allow one to 
split the source galaxies into successively 2 and 5 bins. For the second survey (space-based like), the source galaxy redshift 
distribution has a median redshift Zmedian ~ 1.5, and we assume a good knowledge of photometric redshifts and split the 
sources into 10 tomographic bins. For each redshift bin i, the weighting function is given by 

'^•+1 , , y' 



giix) = 



X 



X' 



0, 



X < Xi+1 

X > Xi + l 



(16) 



where ni{x) is the bin normalized redshift distribution. The average number density of galaxies in this bin is $in with the 
fraction $i given by 

= dx'nix'). (17) 

For example, in the 2 bin case, the normalized distributions are given by 
( n(z) 



(18) 



0, 



Zp > 2zo, 
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1 00 1 000 1 0000 

/ 

Figure 4. The improvement of the constraints as a function of the scales probed by the ground-based hke weak lensing survey (no- 
tomography). The combination PLANCK1-|-SN[0.8]+WL is used here. The curves represent <^ pi^ancki+sn[o.is\/^ plancki+sn[o.is\+wl 
for uiQ and w\. The improvement from WL arises from probing non-linear scales {I > 100), with a significant jump in the improvement 
between I ~ 200 and ~ 2000 for the reference survey considered. 



and 

B, . Jo, Zp < 220, 

For the 5 bins, we use redshift intervals of = 0.6 over the redshift range 0.0 < 2 < 3.0, and for the 10 bins we use intervals 
of Az = 0.3. The convergence power spectra for the 10 bins with the respective sample variance errors averaged over bands 
in i are shown in figure 



3.3 Fisher matrices for weak lensing 

If the convergence field is Gaussian, and the noise is a combination of Gaussian shape and instrument noise with no intrinsic 
correlations, the Fisher matrix is given by: 

where the uncertainty in the observed weak lensing spectrum is given by: l |Kaiser (1992)] [Kaiser (1998)} 

where fsky = 0^7r/1296OO is the fraction of the sky covered by a survey of dimension Q in degrees, and {'yfnt)^^^ is the 
intrinsic ellipticity of galaxies. Guided by previous studies and taking into consideration the major difficulty of constraining 
both Wo and wi, we considered an almost full sky ground-based like survey with f^ky ~ 0.7, a median redshift of roughly 1, 
an average galaxy number density of n = 30 gal/arcmin^, and {'~fi„t)^^^ = 0.4. We also model an ambitious space-based like 
survey with fsky = 0.01, 0.10, and 0.70, a median redshift of roug hly 1.5, n = 100 gal/arcmin^ and {^fntY^^ « 0.25. We 
have used £nis.x = 3000 since on smaller scales, the assumption of a Gaussian shear field underlying Eq. 1201 and the halofit 
approximation to the nonlinear power spectrum may not be valid for larger I's. For the minimum I, we take the fundamental 
mode approximation: 



i... ^ ^ = ,f^, (22) 

f V /"fey 

i.e. we consider only lensing modes for which at least one wavelength can fit inside the survey area. For tomography, the 
Fisher matrix is generalized using 



F^, = + l/2)/.,,Tr fc,-|^C7^gi) , (23) 

withiVf"' =5««'(7L)/*»n 



the power spectrum of the noise in the measurement. 
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PLANCK+SN(z,„„='(5.8l+WLT21t;|.'=0.7J 
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Figure 5. Tomography and the equation of state parameters: the la confidence two-dimensional regions (Ax^ = 1) for the parameters 
{wo,wi). The plots contrast the constraints from different combinations of 1 year of data from Planck, 2000 uniformly distributed 
supernovae with Zmax = 0.8, 1.5, the ground-based like lensing survey (WL), the ground-based like survey with successively 2 bin 
tomography (WLT2) and 5 bin tomography (WLT5), and the deep space-based like survey with 10 bin tomography (WLTIO). 



3.4 Weak lensing systematic effects 

The probes considered here have systematic errors and nuisance factors that need to be well understood and well controlled 
in order for these constraints to be achievable. For cosmic shear, several systematic effects have been identified so far, see 
jRefregier (2003)1 1 ^^'^ references therein for an overview. 

In this analysis, we included the effect of the shear calibratio n bias l|Erben et al.(2001)| [Bac oii et al. (2001)| l 
jHirata fc Seljak (2003)||Bernstein fc Jarvis (2002)1|Kaiser (2000)||Van Waerbeke fc Mellier (2003)t on our results by marginal- 
izing over its parameters. In this bias, the shear is systematically over or under-estimated by a multiplicative factor, and mimics 
an overall rescaling of the shear power spectrum. We used the absolute power calibration parameter ("s and the relative cal- 
ibration parameter ("r between two redshift bins, following the parametrization used and discussed in jlshak et al. (2004)| |. 
The calibration bias is not detected by the usual systematic tests such as the E and B modes decomposition of the shear field 
and the cross-correlation of the shear maps against the point-spread function (PSF) maps. In a weak lensing survey, and 
are parameters of the experiment that can in principle be determined by detailed simulations of the observations. We impose 
in this analysis a reasonable Gaussian prior of 0.04 on these parameters. 

Another serious systematic is the incomplete knowledge of the source redshift distribution | |Wittman et al.(2000)| | 
jlshak fc Hirata (2005)| l including the redshift bias and scattering. One can argue that as long as the scatter in the redshift 
(17(2) « 0.05) is much smaller than the width of the redshift bin Az = 0.3 (for our 10 bin tomography) the effect on the 
integrated results should be small. This is due to the fact that the scatter will only change slightly the shape of the edges 
of the window function that is used for the redshift integration within each bin. At the contrary, the redshift bias alters 
the overall distribution and has been known to affect significantly the cosmological parameter estimation l |Refregier (2003)"] 
[Ishak et al. (2004)| ITereno et al. 20041 . A remedy to this poor knowledge of the redshift distribution using spectroscopic 
redshift has been explored recently in l |Ishak fc Hirata (2005)| |. In the present analysis we marginalize over the redshift bias 
by including the characteristic redshift of the distribution as a systematic parameter Zs and we assume a reasonable Gaussian 
prior of 0.05 on this parameter. 

It is important to question our assumption of gaussianity for the two weak lensing systematic errors that we considered 
here. So far, studies have made this assumption. The motivation for this assumption is the simplicity of the aproximation but 
its justification needs to be addressed in dedicated lensing systematic studies. Thus our treatment of these two systematics is 
only valid under this assumption of gaussianity. For the redshift bias, one has to stress the requirement for a sufficiently large 
number of spectroscopic redshifts, see e.g. l |Ishak fc Hirata (2005)| l, and a sufficiently large number of source galaxies in order 
to reduce the uncertainties to a point where they can be approximated by Gaussians. Also, more narrow band colors and more 
accurate magnitudes (i.e. deeper exposure) are necessary in order to break degeneracies between the photometric redshifts. 
For the shear calibration bias, the errors cannot be made small enough by adding more data, and in this case simulations are 
necessary in order to estimate and address the non-gaussianity question. Massive sky image simulations on which the shear is 
measured and compared to the input will be necessary and such studies are planned to be carried out within ongoing projects 
in the lensing communauty, such as the Shear TEsting Programme (STEP) project jSTEP collaboration (2005)}. Indeed , 
studying the shear calibration bias using massive simulations is one of the goals of this project i js^^^ToUabo^^ion (2005')] l . 
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Hence, the reader should be aware that our treatment of shear cahbration errors is vahd only under the gaussianity assumption, 
and this needs to be checked by future simulation studies. 

Finally, it is important to recall that there are systematic effects that were not included in this analysis (e.g. intrinsic align- 
ments of galaxies, uncertainties associated with the non-linear mapping of the matter power spectrum, see ( [Refregier (2003)] l 
and references therein for a list of other effects), however it is encouraging to note that many efforts are spent in order to 
study these and other lensing systematic effects and some progress has been made. With a better understanding of these 
limiting factors it will be possible to parameterize them and evaluate their effect on the cosmological parameter estimation. 



4 PROBING DARK ENERGY WITH SUPERNOVAE TYPE lA 

Supernovae of type la are powerful probes of dark energy, as when properly calibrated they become cosmological standard 
candles that can be used to measure distances as a function of redshift. The luminosity distance to a SN la is given by 



where L is the intrinsic luminosity and F is the observed flux. The apparent magnitude of this SNa la can be written as 

m = ^\og^a{DL)+M (25) 

where Dl = HqcIl/c is the dimensionless luminosity distance, A4 = — 5 logj^Q{Ho/c) + constant is the magnitude parameter, 
and M is the absolute magnitude, degenerate here with the Ifubble parameter. In a spatially flat model 

DLiz) = {l + z) [ ^ =dz', (26) 

where £{z) is as deflned in Eq.J^J. We use the Fisher matrix for the SNe la defined as (see, e.g. ( |Tegmark et al.(1998) | 
[Huterer fc Turner (200T)| |) 

p - 1 ^^L.i dDL,i , . 

' fri<^n.{Dr.,.r dp'^ dp0 ■ ^^'> 

We use two sets of 2000 SNe la uniformly distributed with Zmax ~ 0.8 and Zmax = 1.5. It is important to briefly note here that 
there are systematic uncertainties associated with supernova searches: these include luminosity evolution, gravitational lensing 
and dust; see, e.g. | |SNAP (2003b)| l and references therein. In order to partly include the effect of these systematics and the 
effect of the supernova peculiar velocity uncertainty l |Tonry et al. (2003)| l, we follow \Kim et al. (2003)||Upadhye et al. (2004)^ 
and use the following expression for the effective magnitude uncertainty 



"■+(1;^) >^?" (28) 

where, a peculiar velocity of (t„ = 500km/sec is assumed, and following ( |Kim et al. (2003)] [Upadhye et al. (2004)| | we use 
5m = 0.02 for space-based supernova survey data, and assume 5m ~ 0.04 for ground-based survey data. The quadrature 
relation I28II assures that there is an uncertainty floor set by the systematic limit 5m so that the overall uncertainty per bin 
cannot be reduced to arbitrarily low values by adding more supernovae. 



5 PROBING DARK ENERGY WITH CMB AND COSMIC COMPLEMENTARITY 

CMB is a powerful cosmological probe; however, like other probes it suffers from some parameter degeneracies and needs to be 
combined with other data sets in order to provide tight constraints on dark energy parameters. For example, it is well known 
that cosmic shear and CMB have different types of degeneracies in their parameters, which are nicely broken when these 
probes are combined. Indeed, among the orthogonal directions of degeneracy between cosmic shear and CMB, are the known 
doublets {Q,m,o-a), {h,ns), and {ns,as); see, e.g. l|Tereno et al. 200*4^ . We use this cosmic complementarity in the present 
analysis where the statistical error (with some systematics included) on a given parameter is given by: 

g\p°') ^ [{¥cMB + FwL + Fsive + n)-^", (29) 

where fl is the prior curvature matrix, and Fcms, Fvkl and Fgjve are the Fisher matrices from CMB, weak lensing, and 
supernovae, respectively. We only impose priors on the characteristic redshift and the calibration parameters by taking priors 
of (j{zp) = 0.05 and (j{C,s) = o"(Cr) = 0.04 on the calibration parameters (corresponding to 2% rms uncertainty on the 
amplitude calibration; jHirata fc Seljak (2003)^ ). We may add the CMB and weak lensing (WL) Fisher matrices together 
because the primary CMB anisotropies are generated at much larger comoving distance than the density fluctuations that 
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Table 2. Summary of the parameter {wo,w-i) estimation errors (Ict uncertainties) from different combinations of probes. CMB exper- 
iments are WMAP 8 years, Planck 1 year, and ACT 1 year (unlensed spectra) combined witli WMAP 8 years. (We did not include 
constraints from the Sunayev-Zeldovich effect or lensing of the CMB that ACT and PLANCK will be useful with.) The 2000 supernovae 
are uniformly distributed with Zmax = 0.8. WL is for a ground-based like weak lensing survey with /^fej, = 0.7, n = 30 gal/arcmin'^, 
i'^int)^^^ ~ 0.4, and a median redshift z^^^ « 1. WLT2 refers to the same weak lensing survey but with 2 bin tomography. The best 
constraints from the combinations in this table are from PLANCK1+SN[0.8]-|-WLT2. 





CMB 


+SN 


-fWL 


-I-WLT2 


+SN-I-WL 


+SN-I-WLT2 


SN+WL 
no-CMB 


SN+WLT2 
no-CMB 


a{wo) 
a{wi) 


WMAP-8 alone 
3.73 
5.65 


0.25 
0.87 


0.66 
1.45 


0.35 
0.66 


0.21 
0.59 


0.11 
0.26 


0.24 
0.93 


0.11 
0.35 


a{wo) 


WMAP8-fACT 
0.82 
1.87 


0.20 
0.59 


0.58 
1.40 


0.22 
0.42 


0.18 
0.48 


0.11 
0.25 






a{wo) 
o-(wi) 


PLANCK- 1 alone 
0.50 
0.31 


0.17 
0.23 


0.28 
0.20 


0.23 
0.18 


0.093 
0.083 


0.086 
0.069 







give rise to weak lensing, hence it is a good approximation to take them to be independent. For CMB, we project constraints 
from 8 years of WMAP data (WMAP-8), 1 year of PLANCK data (PLANCK-1), and 1 year of ACT data combined with 
WMAP-8. The experiment specifications used for Planck and ACT are listed in table For 8 year WMAP, we include TT, 
TE, and EE power spectra, assuming fsky = 0.768 (the KpO mask of Ref. ( [Bennet et al. (2003)} ), temperature noise of 400, 
480, and 580 /iK arcmin in Q, V, and W bands respectively (the rms noise was multiplied by y/2 for polarization), and the 
beam transfer functions of Ref. (|Page et al. (2003)]|. 



6 RESULTS AND DISCUSSION 

We calculated future constraints on dynamical dark energy parameters obtained form several combinations of cosmic microwave 
background experiments (CMB), supernova searches (SNe la), and weak lensing surveys (WL) with and without tomography. 
For CMB, we considered future constraints from 8 years of data from WMAP (WMAP8), one year of data from Planck 
(PLANCKl), and one year of data from the Atacama Cosmology Telescope (ACTl). We used two sets of 2000 supernovae 
with Zmax = 0.8 (SN[0.8]) and Zmax ~ 1-5 (SN[L5]) respectively, and considered various cosmic shear reference surveys: an 
almost full sky (70%) ground-based like survey (WL) with successively 2 (WLT2) and 5 (WLT5) tomographic bins; a deep 
space-based like survey with 10 tomographic bins (WLTIO) covering successively 1%, 10% and 70% of the sky. We combined 
these experiments in doublets and triplets, taking into account space-based or ground-based like experiments for supernovae 
and weak lensing. One should take the uncertainties obtained on the dark energy parameters from CMB-only Fisher matrices 
with some caution as some of them are large and the Fisher matrix approximation may not be valid. We compared our 
CMB-only constraints for Planck with those of reference (| Hu (2001)} and found them in good agreement. For example, when 
we fix the parameters wi to compare with reference (|Hu (2001)} , we find that our {a{^lA) ~ 0.087, a{w) = 0.31} are in good 
agreement with {a{Q,A) = 0.098, a(w) — 0.32} from | |l^^20^^^ Importantly, the constraints we obtain from any combination 
of the three probes are significantly smaller then CMB-only and therefore the constraints obtained are good estimates of the 
low bound of the uncertainties. Our results are summarized in tables ^ O and |1] and in figures |5| |3 |1] and 13 We looked 
for combinations of experiments that will provide constraints that are small enough to answer conclusively some of the dark 
energy questions. 

The first question that one would like to answer is whether dark energy is a cosmological constant or a dynamical 
component. The most decisive answer will be to rule out significantly the cosmological constant. A less decisive but very 
suggestive answer will be to show that dark energy parameters are those of a cosmological constant with a few percent 
uncertainty only. This can be compared to the case of spatial curvature in the universe: when CMB results constrained the 
total energy density to fir ~ 1.02 ±0.02 at the one sigma level, then it became generally more accepted that spatial curvature 
is negligible. However, it is important to note that it will remain always possible to build dark energy models that could have 
a set of parameters indistinguishable from those of a cosmological constant within the limits set. So in this scenario other 
types of tests than the equation of state will be required in order to close the question. 

As shown in tableU the combinations PLANCKl-hSN[1.5]-hWLT10[/ss,a = 0.1] and PLANCKl-hSN[0.8]+WLT5[/sfca = 
0.7] provide impressive constraints that reach the goal set. This is followed by the constraints from PLANCKl-f SN[1.5] + 
WLT10[/sfcy = 0.01]. One should note here that for the equation of state parameters, only small additional improvements to 
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Table 3. Same as table l2l but for the dark energy parameterization {wo,Wa)- As usual the errors on Wa are larger than those on uii 
(roughly equal to twice the errors on uii). The best constraints in this table are from PLANCK1+SN[0.8]+WLT2. 





CMB 


+SN 


+WL 


+WLT2 


+SN+WL 


+SN+WLT2 


SN+WL 
no-CMB 


SN+WLT2 
no-CMB 


ct(wo) 


WMAP-8 alone 
1.84 
3.03 


0.24 
1.25 


0.46 
1.61 


0.39 
1.24 


0.21 
0.92 


0.14 
0.53 


0.26 
1.37 


0.14 
0.76 


(t(wo) 


ACT+WMAP8 
0.52 
1.83 


0.21 
0.92 


0.27 
0.85 


0.42 
1.48 


0.19 
0.76 


0.14 
0.50 




J) 
!) 


t(wo) 


PLANCK-1 alone 
0.67 
0.52 


0.15 
0.30 


0.28 
0.32 


0.24 
0.29 


0.097 
0.133 


0.088 
0.111 




35 



these constraints are obtained when an extremely ambitious sky coverage of 70% is considered for PLANCK1+SN[1.5]+WLT10 
[fsky = 0.7]. Finally, the constraints on the equation of state parameters from PLANCKl+SN[0.8]+WLT2[/sfej, = 0.7] are not 
small enough for the criterion of a few percent set above. 

Another test to answer the same question above is by probing directly the dark energy density at various redshift 
points. For example, if future data will show significant departures of the parameters fi or £2 (see Sec. l2.2L from unity then 
a cosmological constant can be ruled out conclusively. Our analysis shows that this might be a less difficult test as even 
PLANCK+SN[0.8]+WLT5[/,fcy = 0.7] has the potential to achieve ai£i) = 0.012 and (7(^2) = 0.049. 

A second question of interest is what combination of experiments could distinguish between some currently proposed 
models of dark energy. Of particular interest are models that predict an equation of state with parameters that are sig- 
nificantly different from those of a cosmological constant. For example, one could consider quintessence tracker models 
jZlatev et al.(1999)| l or super-gravity inspired models ( |Brax et al.( 19 997), f^'' which, wq > —0.8 and dw/dz{z = 0) ~ 0.3. 
From our tables, we see that the distinction between these models and a cosmological constant can be achieved by several 
different combinations of experiments with different levels of precision. 

Our tables 121 and 13 show that after combining Planck and supernova constraints, weak lensing without tomography adds 
an improvement of roughly a factor of 2 or better to the constraints. As shown in figure 21 the WL-improvement arises 
from probing non- linear scales (/ > 100), with a significant jump between I ~ 200 and ~ 2000 for the ground-based like 
survey considered. Adding 2 bins tomography to the lensing survey provides an additional factor of 2 improvement to the 
combination WMAP8-|-SN[0.8]-|-WL and to combination ACT1-|-WMAP8-|-SN[0.8]-|-WL. We mention here that these results 
do not include constraints from Sunayev-Zeldovich effect or lensing of the CMB, that ACT and PLANCK will be useful with. 

In table 21 and figure |3 we summarize our results on multiple-bin tomography. The constraints on the equation of state 
parameters from PLANCK1-|-SN[1.5] improve by factors 3-5 when WLTIO with fsky ~ 0.1 is added to the combination. 
Also, the constraints on the equation of state parameters from PLANCKl-|-SN[0.8]-|-WLT5[/sfej, = 0.7] are roughly factors of 
3-6 better than those obtained from PLANCK1-|-SN[0.8]. Thus, we find that the improvements obtained from multiple-bin 
tomography lensing surveys are important for the questions raised at the beginning of this section as it brings the uncertainties 
sigiuficantly closer to the goal of a few percent. Therefore, the present study shows that tomography is very useful to add 
further improvements to the constraints on dark energy parameters using both ground-based experiments and space-based 
experiments. The precise discussion of the technical and instrumental feasibility of multiple-bin tomography from ground or 
from space is beyond the scope of this paper and should be addressed somewhere else. 

We took care to include some systematic effects in our analysis. We parameterized the weak lensing calibration bias and 
assumed reasonable priors of 0.04 on the calibration parameters. We also parameterized the redshift bias and used a reasonable 
prior of 0.05. However, there are other systematic effects that we did not include and that may affect our results. For example, 
we did not include the effect of intrinsic correlations between the lensing source galaxies on our results jCroft fc Metzler (2000) | 
[Heavens et al. (2000) | [Lee fc Pen(2'000]1|Catelan et al.(2001)| [Crittend en et al. (2001)] [Brown et al.(20027^ 
jJing (2002)[[Jarvis et al.(2003)[ [Heymans et al.(2004)[[Hirata fc Selja k (2004yf, and we used the HALOFIT fitting formula 
to evaluate the non-linear matter power spectrum (full simulations should be used for real data analysis). As discussed in 
section 3.4, our inclusion of these two lensing systematics assumed their gaussianity. The effect of these systematics on our 
results is thus valid only under this approximation. Future dedicated studies of weak lensing systematic effects should address 
the issue of non-gaussianity. Also, following previous work, we included in our supernova constraints a conservative systematic 
limit, but more studies need to be done there as well. 

Nevertheless, our results show that after the combination of CMB, supernovae, and weak lensing surveys, tomography 
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Table 4. Various tomography analyses: a comparative summary of the constraints (Icr uncertainties) on the different dark energy 
parameterizations from Planck-1, 2000 uniformly distributed supernovae with Zmax = 0.8, 1.5, the ground-based like lensing survey with 
successively 2 bin tomography (WLT2) and 5 bins (WLT5), and the deep space-based like survey with 10 bin tomography (WLTIO). 
The results are presented for the dark energy parameters {wo,w\}, {wo,Wa} and {£\ = pdei^ = 0-5)/pde (0)i ^^2 = pdei^ = l-0)/Pde{0)} 





PLANCK-1 


+SN(z™„^ = 1.5) 


+ SN{z„„^ = 0.8) 


-l-SN(0.S)-|-WLT5 


-l-SN{1.5)-l- WLTIO 


-t-SN{1.5)-l- WLTIO 


+SN(1.5)+WLT10 




alone 




-t-WLT2[/,fc^ = 0.7] 


Ihky = 0.7] 


If sky = 0.01] 


Ifsky = 0.1] 


Ifsky = 0.7] 




0.50 


0.11 


0.086 


0.04 


0.048 


0.032 


0.023 


o-(ii)i) 


0.31 


0.13 


0.069 


0.034 


0.042 


0.027 


0.021 




0.67 


0.12 


0.088 


0.041 


0.049 


0.033 


0.023 


(T(Wa ) 


0.52 


0.21 


0.111 


0.056 


0.067 


0.040 


0.026 


a(£i) 


0.11 


0.048 


0.029 


0.012 


0.013 


0.010 


0.009 




0.32 


0.12 


0.065 


0.049 


0.082 


0.040 


0.018 



with very large fractions of the sky and many redshift bins has the potential to add key improvements to the dark energy 
parameter constraints by bringing them to the level of a few percent. 

On the other hand, the requirement for very ambitious and sophisticated surveys in order to achieve some of these 
constraints, and the difflculty to obtain any further significant improvements, even with the most ambitious survey we 
considered, suggest the need for new tests to probe the nature of dark energy in addition to constraining its equation of state. 
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